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ABSTRACT 

The Be/X-ray binary system A0535+262 underwent a giant outburst in May-June 2005, followed by 
a dimmer outburst in August-September 2005. This increased intensity provided an opportunity to 
search for redshifted neutron-capture lines from the surface of the neutron star. If discovered, such 
lines would constrain the neutron star equation of state, providing the motivation of this search. The 
spectrometer (SPI) on board the INTEGRAL satellite observed the dimmer outburst and provided 
the data for this research. We have not detected a line with enough significance, with the width- 
dependent upper limits on the broadened and redshifted neutron capture line in the range of (2-11) 
x 10 -4 photons cm~ 2 s _1 . To our knowledge, these are the strongest upper limits on the redshifted 
2.2 MeV emission from an accreting neutron star. Our analysis of the transparency of the neutron 
star surface for 2.2 MeV photons shows that photons have a small but finite chance of leaving the 
atmosphere unscattered, which diminishes the possibility of detection. 

Subject headings: nuclear reaction, nucleosynthesis, abundances - stars: neutron - X-rays: binaries - 
stars: individual (A0535+262) 



1. INTRODUCTION 

Matter that is being accreted onto a neutron star has 
such large energies that it can undergo nuclear reactions 
by colliding with particles in the star's atmosphere. Nu- 
clei that are heavier than Hydrogen will therefore break 
apart and generate free neutrons. The nuclear spalla- 
tion p rocess of 4 He, discussed in detail in iBildsten et al.1 
(1993) (hereafter, BSW), be gins with the 4 He atom en- 
tering the star's atmosphere with a very high kinetic en- 
ergy. The atom undergoes many Coulomb collisions with 
atmospheric electrons and slows down. The atom may 
hit an atmospheric proton and splits into a neutron and 
3 He. After this initial spallation, the 3 He travels deeper 
into the atmosphere and thermalizes at a high column 
density. The thermal 3 He can do one of three things: a) 
absorb a neutron (the most possible outcome) , b) collide 
with a fast proton and be destroyed, or c) participate in 
a fusion reaction. The neutron that was liberated during 
the initial spallation undergoes elastic scatterings with 
atmospheric protons, thermalizes and drifts down due to 
gravity. Since most of the 3 He atoms will capture the 
free neutrons, the rate of production of 2.2 MeV photons 
will depend on the amount of excess neutrons. Another 
restraining factor is the depth to which neutrons travel 
before recombining with protons. If recombination (and 
therefore 2.2 MeV photon production) occurs deep in 
the atmosphere, the photon will have to undergo many 
Compton scat t ers b efore it can leave the atmosphere. 
IBildsten et al.l (|1993l ) have shown that for a neutron star 
with moderate magnetic field strength (~ 10 9 G), a ther- 
mal neutron drifts approximately 1 Compton length be- 
fore recombining and have calculated that only 10% of 
the 2.2 MeV photons will escape unscattered. 
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These photons, emitted with an intrinsic energy Eq 
from the atmosphere of a neutron star of radius R and 
mass M, experience a gravitational redshift for an ob- 
server at infinity to an energy E given by 

R _ 2G 

(|Ozel fc Psaltisll2003D . Measuring E/E„ would allow 
a direct measurement of the gravitational redshift and 
place significant constraints on the underlying equation 
of state. Using the P-p relation for a given model, the 
equation of hydrostatic equilibrium can be numerically 
integrated to obtain the global mass-radius (M-R) rela- 
tion for that equation of state, de fining a single curve 
in th e neutron star M-R diagram (|Lattimer fc Prakashl 
l2004h . While the ratio M/R cannot identify the nuclear 
equation of state uniquely, it can rule out several mod- 
els and place constraints on the remaining models. If the 
mass can be measured independently, using other means, 
then the nuclear equation of state can be uniquely iden- 
tified. 

When calculating the 2.2 MeV photon production 
mechanism, BSW have assumed a moderate magnetic 
field and that matter arrives radially onto the neutron 
star. The magnetic field o f A0535+262 has rece ntly been 
measured as ~ 4 x 10 12 G (jCaballero et al.ll2007f) . In such 
high magnetic fields, the ability of such line photons to 
escape the neutron star is contingent upon the optical 
depths for Compton scattering and also for magnetic pair 
creation being less than unity. The Compton scattering 
constraint amounts to one on proximity of the emission 
point to the upper layers of the stellar atmosphere; as 
mentioned above, it was addressed in BSW, using non- 
magnetic Compton scattering cross sections. The pair 
production constraint is more involved, and will be ad- 
dressed at some length in the Discussion below. The bot- 
tom line is that pair creation transparency of the mag- 
netosphere to 2.2 MeV photons will be achieved if their 
origin is at polar magnetic eolatitudes less than around 
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Table 1 

Mass and Radius Values for Common Equations of 
State for Neutron Stars' 1 and the Expected 
Energies of the 2.2 MeV Line 
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1.41 
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1.30 


TI 
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12-15 


1.56-1.72 


MF 


2.7 


13.5 


1.41 



a Shapiro & Tcukolskv (1983) 
b tt: Pion Condensate, R:Reid's, BJ: Bethe-Johnson, TNI: 
3-Nucleon Interaction Approximation, TI: Tensor Interaction, 
MF: Mean Field Approximation 



20 degrees and they are beamed within around 15-20 de- 
grees relative to their local surface magnetic field vector. 
These two requirements limit the expectation for visibil- 
ity of a 2.2 MeV emission line to small, but non-negligible 
probabilities. 

The 2.2 MeV line is expected to be both redshifted 
and broadened, due to gravitational and relativistic ef- 
fects (jOzel fc Psaltisl[200l . Since we don't know a priori 
how much the line will be redshifted, we scanned the en- 
tire 1 MeV - 2.2 MeV region - wide enough to contain 
the redshifted line according to several different equa- 
tions of state ()Shapiro & Teukolskvi ri983) given in Table 
1. The width of the redshifted 2.2 MeV line is also un- 
known, which led us to scan the 1-2 MeV region using 
different energy bands: 5keV, llkeV(AE/E = 0.5%), 22 
keV (AE/Eo = 1%), 44 keV (AE/E = 2%), 100 keV, 
150 keV and 170 keV. We expect the line broadening 
to be minimal, since it is not a rapidly rota ting neutron 
star w ith a pulsation period of 103 seconds (jClark et al.l 
fl998l) . 

The redshifted 2.2 MeV line from A0535+ 262was 
previously studied by iBoggs fc Smith! (|2006h . using 
RHESSI data. The RHESSI data was gathered during 
A0535+262's giant outburst (May 18 - June 25), whereas 
we used data from a smaller ou tburst that took place ap- 
proximately 3 months later (see lCaballero et"aT1l2007l for 
a detailed discussion of INTEGRAL and RXTE analysis 
of the source). The RHESSI instrument uses Germa- 
nium crystals as detectors like SPI, but its area is smaller, 
and there is no shielding. It is specifically designed to 
study solar flares, and due to the close proximity of the 
Sun to this source during the giant outburst, while IN- 
TEGRAL wasn't able to study this event, RHESSI had 
an unique opportunit y to observe the entire outburst. 
IBoggs fc Smith! (|2006D subsequently calculated flux up- 
per limits for A0535+262 which we used for comparison. 

2. SPI INSTRUMENT AND ANALYSIS 

SPI is the spectrometer on board the INTEGRAL 
satellite that was launched in October 2002. SPI con- 
sists of 19 high purity Ge detectors, although only 17 are 
currently operational since the loss of detectors 2 and 17 
on December 6, 2003 and July 18, 2004 respectively. SPI 
has an energy range of 20 keV - 8 MeV with an energ y 
resolution of 2.3 keV at 1.3 MeV (jVedrenne et al.ll2003r ). 
This high precision makes SPI a good candidate for de- 
tecting narrow lines. The spectrometer uses coded aper- 
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Fig. 1.— The light curve of A0535+262 for the interval 
30.01.2005 - 26.11.2005. The spike around MJD 53530 represents 
the giant outburst and the smaller feature around MJD 53620 rep- 
resents the normal outburst in August - September 2005. 



ture for imaging and is shielded by an anti-coincidence 
system that issues a veto flag for every photon that hits 
the detectors without passing through the mask. 

The INTEGRAL data we used was gathered between 
August 31 and September 2 2005 (revolution 352). The 
good observing time was about 209 ksec, with a to- 
tal of 56 pointings. We tried four different background 
subtraction methods: two different empty fields, GED- 
SAT and GEDRATE. The first empty field we used was 
from revolution 220 (standard empty field in OSA, ob- 
servations of Algol) and consisted of 103 pointings and 
186 ksec. The second empty field we used was from 
revolutions 348 and 349 (observations of Cen-A, cor- 
responding to August 2005) closer by date to our ob- 
servation, and consisted of 41 pointings and 123 ksec. 
Although the latter one wasn't really "an empty field" 
in hard X-rays, no source was detected above 1 MeV. 
The GEDSAT and GEDRATE background models are 
integrated in the SPI analysis software. GEDSAT as- 
sumes that the background level is proportional to the 
product of saturated detector trigger rate and live-time, 
whereas the GEDRATE background model assumes that 
the background level is proportional to the product of 
non-saturated detector trigger rate and live-time. 

We used spiros version 5.1 for data analysis, provided 
by ISDC. For each of the four background subtraction 
methods we accumulated data for detector subsets: sin- 
gle events (SE), PSD events (PE), and double events 
(ME2). The single events represent the cases where a 
photon is detected by any one of the 17 are currently op- 
erational, but not tagged by the PSD. The PSD events 
represent the events which are analyzed by the Pulse 
Shape Discriminator. The double events (or multiple 
events with two detectors, ME2) are defined as cases 
where groups of two detectors detect an event simul- 
taneously, in which case the energy of the photon is 
given by the sum of the two detectors. The reason we 
have different detector subsets is as follows. We do not 
want to use only SE+PE events (which is usually done 
for analysis of low energy data) because there is sig- 
nificant area at ME2 above 1 MeV. Moreover, for SE 
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there is a si gnificant background noise in the 1.4 - 1.7 
MeV region (jWeidenspointner et al.l [2003) which is arti- 
ficial. These background features are diminished when 
double events and/or PSD events are used. Therefore 
except in this region we merged PE, SE, and ME2, while 
within the region of the artificial background component 
we used PE+ME2. This was done by using a modified 
spi-obs-hist, and renormalizing the effective area after 
the fluxes are obtained. The rest of the analysis has 
been done using standard OSA software. 

In order to create the spectrum of A0535+262 we 
did imaging in 5keV steps and plotted the flux (pho- 
tons cm -2 sec -1 keV -1 ) against energy. We used IDL 
to fit Gaussian curves of given energy and FWHM and 
extracted the errors of the flux, given by the area under 
the Gaussian, to calculate the 3 a flux upper limits. 

3. RESULTS 

Two spectra of A0535+262 are given in Figures 2 and 
3, for different background subtraction methods and dif- 
ferent detector subsets. The significance of A0535+262 
did not exceed a few a in any of the imaging analyses, 
therefore we couldn't detect a line. However, we calcu- 
lated upper limits for the flux for different background 
subtraction methods and line widths. In order to do 
this, we fit the residual spectrum to Gaussian lines with 
a constant background in 5 keV steps. The errors of 
these fits provided us with the flux upper limits. When 
fitting Gaussians with large FWHM values (such as 44 
keV or 110 keV), the low and high energy limits of the 
spectrum gave unusually high error values due to the 
lack of data below 1 MeV and above 2.2 MeV; these 
w ere neglected. W e chose to use the same FWHM values 
as iBoggs fc Smith! f|2006T ) in order to make a meaningful 
comparison of the RHESSI and SPI instruments. The 
line widths range from 11 keV FWHM (AE/E = 0.5%) 
to 110 keV FWHM (AE/E = 5.0%) since we cannot 
calculate the exact line broadening, which depends on 
the geometry of the mass accretion and emission. Three 
energy dependent flux upper limit plots are shown in Fig- 
ures 4, 5 and 6 and the upper limits for all background 
subtraction methods and line widths are given in Table 
2 along with the RHE SSI upper limits for comparison 
(|Boggs fc Smithll200l . 

4. DISCUSSION 

Despite the excellent energy resolution and relatively 
high area of SPI for high energy photons, the redshiftcd 
2.2 MeV line from A0535+262 was not detected signif- 
icantly. However, we were able to obtain stronger con- 
straints in the flux than that of RHESSI instrument. 
The non detection is not surprising, as the transparency 
of the neutron star, detailed below, is greatly diminished 
for high magnetic fields, reducing the yield of 2.2 MeV 
photons. 

The visibility of any atmospheric 2.2 MeV neutron 
capture line is subject to potential attenuation in the 
neutron star magnetosphere. The leading process for 
such absorption is magnetic pair production, 7 — > e + e~, 
the process that underpins standard radio an d pola r 
cap gamma-ray pulsa r mod els (e.g. ISturrockl (|1971f ). 
iDaughertv fc Harding) (|1996l )). This mechanism is the 
highest order photon absorption interaction in strong- 
field QED, and is prolific at the field strengths inferred 
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Fig. 2. — The upper panel shows the background subtracted 
spectrum of A0535+262 using revolution 220 as flat field. Only 
single events and PSD events were used in this analysis. The lower 
panel shows the significance of A0535+262 in the 1-2.2 MeV region. 
The 3<t threshold is rarely crossed (not verified by other background 
methods and energy bands), which prompted us to calculate upper 
limits for the flux. 
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Fig. 3. — The upper panel shows the background subtracted 
spectrum of A0535+262 corrected for the SPI effective area and 
using GEDSAT for background subtraction. Single events, double 
events and PSD events were used in this analysis. The lower panel 
shows the significance of A0535+262 in the 1-2.2 MeV region. 



for A0535+262. It has a threshold (for the || polariza- 
tion state) of 2m e c 2 / s'm9kB, where Oub is the angle the 
photon momentum vector makes to the local field line, 
above which the interaction r ate precipitou sly rises (e.g. 
IDaughertv fc Hardind(|l983T ). lBarind(|2008l ). At energies 
only modestly above threshold, the attenuation lengths 
in local fields greater than around 10 12 Gauss are con- 
siderably less than 1cm. Accordingly, once the threshold 
is crossed, 2.2 MeV photons can be expected to be de- 
stroyed in A0535+262, provided this occurs at altitudes 
less than around one stellar radius above the surface. 
This accessibility of the threshold defines the broad cri- 
terion for transparency versus opacity of the magneto- 
sphere for the line photons, the outcome of which is crit- 
ically dependent on geometry. Attenuation phase space 
can readily be assessed as functions of the surface po- 
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Table 2 

Redshifted 2.2 MeV Line, 3 a Upper Limits, 1-2.2 MeV 



AE/Eo (FWHM) 


RHESSI Flux" 


Rev. 220 Flux 6 




Rev 348/349 ft Flux 


GEDSAT Flux 6 


GEDRATE Flux 6 


GEDSAT Flux c 


(%) 


(photons cm -2 s — 1 ) 


(photons cm -2 s~ 


l ) 


(photons cm -2 s^ 1 ) 


(photons cm -2 s — L ) 


(photons cm -2 s — 1 ) 


(photons cm -2 s — 1 ) 


0.5 


4.0 X 10" 4 
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2-3 X 10~ 4 


2-3 x 10" 4 


2-3 X 10" 4 


1.0 


5.0 X 10" 4 
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3-5 x 10" 4 


3-5 x 10" 4 


3-4 x 10~ 4 
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6.4 x 10" 4 
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4-6 x 10~ 4 
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10.5 x 10" 4 
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7-10 x 10~ 4 
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8-11 x 10" 4 
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a Boggs k. Smith ( 2006) 
'Using single and PSD events only. 
c Using single, double and PSD events. 
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Fig. 4. — The 3 a upper limits of A0535+262 using single and 
PSD events, Revolution 220 as flat field and fitting a Gaussian with 
llkeV FWHM (AE/Eo = 0.5%). The peak around 1100 keV is due 
to electron pair production line. The 1400-1700 keV region is domi- 
nated by instrumental background lines, therefore only PSD events 
were taken into account in that region. Th e peak around 1800 keV 
is an other instrumental background line ( Weidcnspointncr ct al. 
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Fig. 5. — The 3 a upper limits of A0535+262 using single, 
multiple and PSD events, GEDSAT for background subtraction 
and fitting a Gaussian with 22keV FWHM. 



lar field, the surface emission locale, and the emission 
direction of the photons. Such a detailed analysis was 
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Fig. 6. — The 3 a upper limits of A0535+262 using single and 
PSD events, GEDSAT as flat field and fitting a Gaussian with 
llkeV FWHM. 



offered by iBaring fc Harding] (|2001[) at much higher en- 
ergies that are germane to pulsar and magnetar contin- 
uum emission. It should be remarked that familiar two 
photon pair creation is negligible in comparison. 

The transparency of the magnetosphere for 2.2 MeV 
line photons was investigated here principally via the 
salient threshold criterion, to quickly ascertain the per- 
tinent approximate phase space. Our analysis was per- 
formed in flat spacetime, sufficient for a general indica- 
tion of transparency, for which rather simple geometri- 
cal relations for photon transport through the magneto- 
sphe re can be prescribed, alon g the lines of the approach 
that IBaring &: Harding] (|2007l ) used for resonant Comp- 
ton upscattering in magnetars. The involved mathemat- 
ical details of this analysis will be omitted here, for the 
sake of brevity, and only the principal conclusions out- 
lined. Photons at 2.2 MeV that are emitted at the atmo- 
spheric surface are below magnetic pair creation thresh- 
old when they initially have sinflfcs2m e c 2 / (2.2MeV) 
OkB < 27°. Note that at the field strengths of < 4 x 10 12 
Gauss appropriate for A0535+262, the pair threshold 
for the _L photon polarization state is extremely close 
in energy to that for the || state. For equatorial emis- 
sion locales, where the radius of field curvature is com- 
paratively small, even photons that are initially beamed 
close to field lines rapidly acquire sufficient angles with 
respect to B to precipitate pair attenuation. Hence the 
magnetosphere is necessarily opaque to 7 — > e + e~ for 
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emission colatitudes 9 e generally greater than around 
30°. Near the pole, if 9kB < 27° initially, transport 
of photons to significant altitudes is required to access 
pair threshold. There, the field line curvature radius is 
p c ks 4R ns /(3 sin# e ) for a neutron star radius of R ns . 
Accordingly, photons that initially are beamed close to 
the field, acquire an angle of 9kB ^ 27° when reaching an 
altitude h ~ AR ns sin 27°/ (3 sin# e ). This altitude must 
generally be more than around a stellar radius for the 
field decline to render the pair creation rate insufficient 
to sustai n photon attenuation, as o pposed to free escape 
(see, e.g. iBaring fc Hardingl (|2001f ) for a detailed exam- 
ination). To summarize our analysis, here we find that 
pair creation transparency of the magnetosphere to 2.2 
MeV photons will occur if their origin is at polar mag- 
netic colatitudes 6 e < 20° and they are initially beamed 
within around 15-20 degrees relative to their local sur- 
face magnetic field vector. These two requirements limit 
the expectation for visibility of any 2.2 MeV emission 
line to modest or small, but non-negligible, surface areas 
and solid angles. 
With a more powerful tool for background rejection 
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and a dramatical improvement in sensitivity (broad line 
sensitivity of 1.2 x 10 -6 ph otons cm" 2 s" 1 and a spec- 
tral resolution of 0.2 - 1 %) (|Boggsll2006l ) . the Advanced 
Compton Telescope could be a good opportunity to fol- 
low up on the search for the broadened and redshiftcd 
2.2 MeV line. 
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